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text. The orders we considered are a superconducting order, modeled by a charged 

scalar field, and a magnetic order modeled by a neutral scalar field. We also discuss 

t — . ■ the case of two competing scalars coupled to a single gauge field. 

In all cases discussed here the phases tend to compete, rather than enhance each 

other. The condensation of one scalar hinders any further instabilities, unless we have 

a sufficiently strong repulsive interactions between the bulk scalars. We provide both 

analytic arguments and numerical demonstration of this fact. 

r> \ Based on the cases discussed here, we conjecture that holographic orders tend to 

compete for attractive bulk interactions, including gravity, and to cooperate, or be 

mutually enhancing, for repulsive bulk interactions between the corresponding order 

parameters. 
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1 Introduction and Conclusions 

Holographic Condensed Matter (for a review see e.g pQ) is a new approach to modelling 
strongly interacting systems and calculating their properties. The method involves a 
set of tools acquired through investigations of quantum gravity, including string the- 
ory, black hole physics and higher dimensional gravitational theories. As such, better 
understanding of the scope and usefulness of the approach for various physical systems 
may improve our understanding of quantum gravity. However, the main utility of this 
approach for the time being is the construction of new classes of strongly interact- 
ing universality classes, which may well be useful to understanding and modeling real 
materials, in cases where conventional calculational tools fail. 

With this phenomenological approach in mind, we discuss in this paper the issue of 
competing orders. The holographic approach, by and large, has concentrated so far on 
the dynamics of a single order parameter coupled to all relevant conserved quantities 
(somewhat similar to the Landau-Lifshitz theory of phase transitions, however not 
necessarily near the phase transition point). Strongly correlated electron systems tend 
to have complicated phase diagrams, with many possible orders, such as magnetic 
orders, striped phases or superconductivity. This is in contrast to conventional Fermi 
liquids, which tends to have very few instabilities. It is interesting therefore to look 
at the interplay of various orders in the holographic context, to model specific phase 
diagrams, and even perhaps to draw some general conclusions. 

In weakly interacting Fermi liquids, the onset of one type of order tends to produce 
a mass gap (at least for parts of the Fermi surface), thus inhibiting further instabilities. 
In this sense the phases tend to compete, the parameter range of a potential instability 
tends to shrink when the competing order sets in, for example the critical temperature 
is lowered. In some cases, the onset of one type of order prevents any further instability. 
This is one reason for the robustness and usefulness of the Fermi liquid picture. 

On the other hand, quantum critical points may change that picture. Such phase 
transitions deviate from the Landau paradigm of phase transitions in that the critical 
theory may have new light degrees of freedom, beyond the order parameter and the 
conserved quantities. These emergent degrees of freedom (for example the emergent 
gauge fields in [2]) can enhance some instabilities. For example, it is speculated that 
new massless modes at a quantum critical point may replace the phonons of the BCS 
theory as the pairing mechanism for superconductivity. 

In this note we initiate a discussion of competing orders in the holographic context, 
and examine a few models. By and large we find that the conventional picture is still 
valid for holographic models. That is, for generic models the onset of one order tends 
to suppress all other potential orders. We do find some exceptions to this rule, and 
comment on the set of conditions that seem to be needed to realize holographically the 
compelling idea of mutually enhancing orders. 

The set of quantum phase transitions discussed here is subject to the general con- 
siderations of [3], in that they are triggered by a bulk field falling below the BF bound, 
and thus are expected to exhibit BKT-like scaling near the phase transition point. 
Some holographic models, with this type of mechanism for instability, indeed show the 
expected scaling behavior |3HS]. It would be interesting to show such behavior in the 
present set of models as well. 

The plan of the paper is as follows: we present the set of models we use in this 
paper in section 2, including a brief review of the realization of superconductivity and 



magnetic order in the holographic context. Section 3 is devoted to the interplay of the 
superconducting and magnetic orders, as realized in the holographic context. Finally, 
in section 4 we investigate similar issues in a model of a single order (global U(l) 
symmetry) coupled to two competing scalar fields. 

2 Background and Setup 

2.1 Holographic Superconductors at Zero Temperature 

Holographic superconductors are gravity backgrounds coupled to a gauge field (repre- 
senting global or weakly gauged U(l) in the dual theory), and a charged scalar field 
whose condensation triggers superconductivity. The specific set of models we consider 
here was constructed by Gubser and Nellore [7] and also by Horowitz and Roberts [8] . 
These models correspond to a 2+1 dimensional CFTs, at finite density and zero tem- 
perature. The presence of chemical potential breaks the conformal symmetry, and the 
theory flows to an IR theory with Lifshitz scaling [9], with dynamical exponent z. The 
limiting case of z = 1 corresponds to relativistic CFT in the IR, whereas the case of 
z = oo is similar to the one discussed in [5] (which discussed the charged black hole 
without scalar fields), in the sense of having an AdS2 space in the IR. 

Consider then the bulk theory described by gravity coupled to a Maxwell field with 
the Lagrangian: 
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where F^ v is the electromagnetic field strength, and Aq = <3?(r) is the electric potential 
(satisfying A(oo) = n at finite chemical potential [i, note that at T = all non-zero 
chemical potentials are related by scaling). The charged scalar ip of charge q triggers 
superconductivity when condensing. When discussing this set of models, we work 
exclusively with the quadratic potential V(ip,ip*) = m^\tp\ 2 with rnL > 0. 
The Lagrangian allows for solutions with Lifshitz scaling, with the metric 

ds 2 = -g( r ) 2 dt 2 + 4(<te 2 + dy 2 ) + e 2B{T) %dr 2 (2.2) 

with 
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*(r) = p~-( 1 -) Z (2.3) 
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where the parameters of the solution are given by 
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Following [7] we parametrize the solutions by a pair [z > 1, nT-3, > 0). There exists 
a solution for each such pair, though such solution is not necessarily obtainable as the 
IR limit of a UV CFT at finite density. Under conditions discussed in detail in [7], 
there are irrelevant perturbations that can be potentially used to connect the solutions 
for appropriate (z,m 2 ) to a complete flow from a relativistic 2+1 dimensional CFT. 
We have obtained numerically such flows for the cases discussed below. 

For later use, let us consider the issue of the Breitenlohner-Freedman (BF) bound 
PEHH] in this set of backgrounds. A neutral scalar field £ satisfies the equation of 
motion 

r 2 i"{r) + (3 + z)r£{r) - L 2 U'(£(r)) = (2.5) 



Therefore the scaling dimension Air for small fluctuations of the field £ satisfies the 
equation 
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A(A-l) + (3 + z)A-L 2 m 2 eff = (2.6) 



where m 2 ff is the curvature of the potential U(£) around the background value of £. 
To avoid complex scaling dimensions, which would signal instability, we have the BF 
bound 

Plugging in the value of Lq from above gives 

■? 3 z + 2 . , 

Note that this expression interpolates between the bound for the case the IR geometry 
is AdSi (z = 1 and therefore m 2 rr > —-^i) and the bound for the case the IR geometry 
is A1IS2 (z = oo and therefore mL > — t^Jj. 

Another set of models we use was constructed by Horowitz and Roberts [8]. In 
these models rrfi = 0, and for large enough charge, there is a superconducting order at 
zero temperature. The geometry in the IR is AdS^ with the same radius of curvature 
as the UV space. The BF bound for a neutral scalar field £ is therefore identical in 
both limits of the geometry. 

2.2 Magnetic order in Reissner-Nordstrom Backgrounds 

In [5], the authors discussed magnetic instabilities in the holographic context. In many 
models of condensed matter physics, the spin rotations decouple from spatial rotations 
at long distances, and can therefore be considered to be a separate global symmetry 
when discussing only IR physics. In the holographic context, such global symmetry 
is modeled by an SU(2) gauge field (which we denote as AV). In order to model 
magnetic order, we are interested in configurations that break this global symmetry 
spontaneously. In particular, an anti-ferromagnetic order does not involve macroscopic 
background spin density, and can be modelled holographically by a scalar field in the 
adjoint ip a . By a suitable rotation, the symmetry breaking can be chosen to point in the 



6 It may seem curious we reproduce the bound for the extremal charged black hole, which has AdS2 in 
the interior but no scalar hair. This is the case because the BF bound depends only on the IR geometry, 
and not on how it is supported. 



3 direction, and we denote the resulting neutral (with respect to the electromagnetic 
U(l)) scalar field i/j 3 = £. The Action for the new fields A a and £ is 
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where 
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the potential for the scalar field £. To avoid clutter, we introduce the mass parameter 
v = mJL 2 . 

The parameter A multiplying the action for the field £ controls the backreaction 
on the geometry. We will work in the probe limit, of A being infinite, thus avoiding 
such backreaction. The action for the excitations of the field £ will be bounded in 
all the backgrounds we consider, thus the probe approximation is justified (unlike an 
analogous approximation for the charged scalar fields). 

In order to introduce finite temperature and density, the authors of |5j work in a 
charged black hole background. This background is described by the metric: 

ds 2 = LV (-/(r)d( 2 + dx 2 + dy 2 ) + £j£- (2.10) 

/W-i + ^-i^ (2.ii) 

and the U(l) gauge potential is Aq = fi(l — ). The temperature and density are related 
to the parameter r\ as \i = y/3rj and T = ^(1 —rj). The extremal limit corresponds to 
rj = 1. 

The extremal Reissner-Nordstrom background is a domain wall, interpolating be- 
tween AdSi in the UV, and an AdS2 space. Since the BF bound is different in the IR, 
there is a window of mass parameter, —2.25 < v < —1.5 where the scalar £ induces 
instabilities in the IR, without destabilizing the complete geometry. In this range, we 
expect there to be a normalizable mode for the scalar £, signaling a symmetry break- 
ing and an anti-ferromagnetic ordering. Such expectation is borne out by numerical 
construction of such normalizable mode, for mass parameters within the interesting 
range. The resulting phase diagram (as function of \x and temperature) is given in [5], 
and we reproduce it below in the appropriate limit of our discussion. 

3 Magnetic Instabilities in Holographic Super- 
conductors 

3.1 Phase Structure without Direct Coupling 

In this section we revisit the issue of magnetic instabilities, in the cases where the IR of 
the theory experiences condensation of charged scalar field, thus giving rise to super- 
conductivity at low enough temperature (and in particular at T=0). We concentrate 
on the case where the IR theory developed dynamical scaling with scaling exponent z. 
The case z = 1 gives rise to 2+1 dimensional conformal field theory in the IR. For the 
case z = oo has dynamical scaling in the time direction only, and the IR geometry is 
identical to the one discussed in [5]. 



In figure ([I]) we draw the BF bound for the neutral scalar £ in the IR geometry. 
We see that the bound depends only on the dynamical exponent z. For any value of 
z > 1 , there is a range of masses for which the scalar £ is unstable in the IR region only, 
and as z increases this window of instability increases. The expectation is therefore that 
magnetic order, in the sense defined in [5], is possible only in the region shown in figure 
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Figure 1: Expected phase diagram, explicit solutions representing holographic supercon- 
ductors were constructed numerically for low values of z. The expected phase diagram is 
confirmed for all those solutions. 



This expectation is borne out by explicit numerical calculations. First, the construc- 
tion of Lifshitz solutions in |7j is an holographic superconductor only when embedded 
in an appropriate UV completion, representing a complete flow. We were able to con- 
struct such flows for low values of z only. At asymptotically large z, the geometry 
becomes close to the extremal charged black hole background, which does not have 
scalar hair. We expect then that for sufficiently large z there are no solutions which 
correspond to asymptotically AdS± flowing to the corresponding Lifshitz interior. In 
figure ([2]) we show the results obtained for the case z = 2. 

For that case (and several other low z cases), we were able to construct normalizable 
solutions for the scalar £, if and only if the mass parameter v is in the instability window. 
For this range of masses, the fluctuations of the scalar £ around the maximum of the 
potential £/(£) becomes unstable in the IR region, and it is therefore close to the 
minimum of the potential, £ = £o- On the other hand, the scalar is stable at the 
UV, therefore it is close to the maximum at £ = 0. The scalar profile is therefore an 
interpolation between these two values, facilitated by perturbing the IR value by an 
irrelevant perturbation, and the UV value by a relevant one. The ansatz for the scalar 
field is therefore 

£{r) = £o + ar AlR (3.12) 



for small r, and 



f( r ) = br Auv 



(3.13) 



where A/# > is the IR dimension of £ around the minimum of the potential, and 
Ajjy < is the UV dimension of the field £ around the maximum of its potential [J 

A few such scalar profiles are displayed in figure ([2]). For generic mass parameters 
within the instability window, the scalar profile lies entirely within the Lifshitz part 
of the geometry. Only when the mass parameter v approaches the lower boundary 
of the instability window, v = —2.25, the profile becomes wider than the Lifshitz 
region, signaling instability of the entire spacetime. On the other boundary of the 
instability window, the solution becomes infinitely thin, and thus within the errors of 
the numerical calculations. 

The case z = 1 corresponds to an emergent conformal symmetry in the IR. In 
such cases the instability window disappears, and we do not expect any normalizable 
solution. Indeed, numerical calculations fail to produce such solution, both for the set 
of background considered by Gubser and Nellore [J] , and those considered by Horowits 
and Roberts [8]. 

The picture emerging from these numerical studies is that once the charged scalar 
condenses and superconductivity kicks in, the bulk geometry changes in such a way as 
to inhibit further instabilities, and in particular it inhibits the instability considered 
by 0. 
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Figure 2: Neutral scalar profiles for z=2 and various values of v within the instability window 
—2.25 < v < —2. The blue curve corresponds to v — —2.23, the red one for v = —2.1, and 
the yellow one for v = —2.01. 



7 For the range of masses considered here, there is a choice of quantization in the UV, and our choice 
corresponds to the one considered in [5] , it would be interesting to consider the alternate quantization, whose 
details a slightly more complex. 



3.2 Adding Direct Couplings 

While the results obtained are fairly restricted, when the two scalars interact only 
gravitationally, much more general results can be obtained when the two bulk scalar 
fields interact directly. Consider for example the modified action for the scalar field £ 
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where rj is a new coupling which can be positive or negative. In the probe limit A — > oo 
this coupling has vanishing effect on the charged scalar tp and the superconductivity it 
expresses. On the other hand, in the IR region (where ip = ipo), this coupling translates 
into a shift of 2t]i/jq in the effective mass of the neutral scalar field £. The condition for 
instability of that scalar in the IR becomes 
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whereas in the UV if) = and thus we still require v > —2.25. It is therefore clear that 
for r\ < we can obtain the situation where superconductivity increases the range of 
masses for which magnetic ordering occurs. Such situation is depicted in figure ([3]). 
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Figure 3: Phase diagram with scalar interaction, and a negative coupling constant (77 
-0.3). 



The situation is complementary for positive values of the coupling 77. For general 
values of z the coupling 77 serves to further inhibit the magnetic ordering represented by 
£, completely diminishing its window of instability for a range of scaling parameters z. 
In particular, at sufficiently low values of z (which is the range more likely to correspond 
to complete superconducting flows) , the direct coupling dominates the effect due to the 
modified geometry, and prevent further instabilities for the scalar £. The phase diagram 
is shown in figure(|4|). 
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Figure 4: Phase diagram with scalar interaction, and a positive coupling constant (rj = 0.2) 



It is curious that in models with emergent CFT [z = 1), such as the ones constructed 
by Horowitz and Roberts and used here, the addition of the quartic coupling discussed 
here, or any other similar couplings, has no effect. In such models superconductivity 
prevents any further ordering from occurring even with the inclusion of direct couplings 
between the order parameters. 

The results in this subsection also clarify the situation occurring in the absence of 
any direct coupling. When the interactions between the holographic order parameters 
are attractive (as in the case of gravity, or a positive direct coupling), we find com- 
petition between the corresponding phases. Only when the interactions between the 
bulk fields are repulsive (such as the interactions mediated by vector fields), can both 
phases can coexist, and even enhance each other. 

We conjecture that this correspondence holds in more general situations. Such 
correspondence could be useful as guide to modeling interesting phenomena for strongly 
coupled electrons, where cooperation between phases is expected to be an important 
part of the story. 



3.3 Superconductivity in Magnetically Ordered Phases 

We now turn to the complementary discussion, how is superconductivity affected once 
the holographic material enters the anti-ferromagnetic phase? 

Suppose we include the backreaction of the scalar £ on the geometry, when the 
coupling A in (|2.9p is now finite. The detailed picture necessitates numerical solution 
of the coupled equations, but the leading order effect of the backreaction on potential 
superconductivity can be estimated by the following argument: 

Since the scalar £ interpolates between its maximum £ = in the UV, and it's 
minimum at £ = £0 m the IR, the leading order change in the geometry would be a 
shift of the cosmological constant in the IR, to be slightly larger than it was without the 
backreaction. The radius of curvature L^ for the resulting AdSi would then decrease 



slightly. 

Consider now the addition of a charged scalar field t/j to the modified geometry, to 
model potential superconductivity. A necessary condition for instability of the ip = 
configuration is [2E] 

- 2.25 < m\L 2 and {m\ - 2q 2 ) h\ < -1.5 (3.16) 

The first condition is the BF bound in the UV, which is satisfied, and the second 
condition is the BF bound in the IR (using the effective IR mass m 2 b — 2q 2 ), which has to 
be violated for the ip = configuration to be unstable, and thus for superconductivity 
to occur. 

We find therefore that the leading order effect of the scalar £ is to modify the 
geometry is such a way as to make the second condition more restrictive, without 
modifying the first one. In other words, for any given charge q, the range of charged 
scalar masses which give rise to superconductivity shrinks due to the presence of the 
neutral scalar £. 

This argument leads to results that are consistent with our general conjecture: for 
gravitational coupling only, or when including direct attractive interactions, the mag- 
netic ordering seems to inhibit superconductivity. However, the addition of repulsive 
coupling may give rise to cooperation, and to superconductivity in the magnetically 
ordered phase, in cases where this would not occur without that magnetic ordering. 
It would be interesting to look at this issue in more detail, either by constructing the 
back-reacted solution, or perhaps by adding direct coupling between the two scalar 
fields. For example, models in which the anti-ferromagnetic phase enhances the pos- 
sibility of superconductivity may be of relevance to the study of superconductivity in 
the cuprates. 

4 Global U(l) with two scalar fields 

In this section we generalize the discussion of holographic superconductors, to include 
two scalar fields coupled to a single Abelian gauge fielqj. The two scalar fields are in 
competition, very much like the examples discussed above, and we will demonstrate, 
using both numerical and an analytic arguments, novel mechanisms by which a con- 
densation of one scalar hinders the condensation of the other. 

We therefore analyze the phase diagram of two scalar fields charged under the 
same U{1) gauge field. This corresponds to having two potential order parameters 
for a global symmetry in the boundary theory. We will ignore their backreaction to 
the geometry since we are primarily interested in the onset of the condensations (i.e. 
we are working close to the point of the phase transition, where the amplitude of the 
condensate is small). The scalars will carry different charges and masses; since they do 
not back-react on the geometry, they are coupled only through the gauge field. 

Consider then the planar limit of the four dimensional uncharged AdS black hole: 

dr 2 1 

ds 2 = -f(r)dt 2 + -P— + Aj(dx 2 + dy 2 ) (4.17) 



r A f{r) r 



^Models of quantum phase transition with two scalar fields were considered in [6]. 



10 



with 

/( r ) = X^2- MT - ( 4 - 18 ) 

The horizon is at r = 1, while the conformal boundary lives at r = 0, L is the radius 
of the anti-de Sitter space and the temperature of the black hole is given by |j 

nn 3M 1 / 3 

In this note we will adopt the convention that M = L = 1. 
The matter field Lagrangian is given by: 

L= f dx 4 ^^(-^F ab F ab + mj^ + m 2 2 ^--\d^ 1 -ie 1 Ai; 1 \ 2 -\d^ 2 -ie 2 Ai; 2 \ 2 ), 

(4.20) 
where A is the U(l) gauge field, and ipi,ip2 are two charged scalar fielda 10 !. 

We will also assume without loss of generality that m\ > in 2 . We will be looking 
for static solutions and will assume all the fields are homogeneous in the field theory 
directions with only radial dependence. 

The equations of motion for the fields in this coordinate system are: 

/' ,. 1 (e 2 A 2 m 2 \ 



/' 1 (e\A 2 m\ 



V>2 + -V>2 + -4 —T2 T ) ^ 2 






A'l - 2e 2 i±A> - 24-p.At = 0. (4.21) 



By rescaling of At we may always set the charge of one scalar to unity and in most 
cases we set e\ = 1. Phase structure of the theory is determined by dimensionless 
ration |^- To require regularity at the horizon we will have to set At = at r = 1. 
Since we have a set of coupled equations, this will in turn give the constraints at the 
horizon 

V4 = ^i 

V>2 = ^2 

A t = (4.22) 

at t = 1. Examining the behavior of the fields near the boundary, we find 

V>1 ~ ^ + A T X +> 1 +'I'_iT A -> 1 + ... 
lp 2 ~ ^ + i2 T A +> 2 +1'_,2T A -> 2 + ... 

A t ~ n + pT d - 2 + ... (4.23) 



9 It should be noted that "z" is used as a co-ordinate here. In the previous section it was used as Lifshitz 
exponent. 

10 From now on, any quantity denned with suffix 1 will be for scalar field ipi and any quantity defined with 
suffix 2 would be for tp2 ■ A quantity related to scalar field defined without a suffix will be a generic quantity 
relevant for both fields. 

11 



The coefficients above can be related to physical quantities in the boundary field theory 
using the usual dictionary in gauge/gravity correspondence. The constants /i, p are 
the chemical potential and the density of the charge carrier in the dual field theory, 

respectively, and X± t i t 2 = \{d ± Jd 2 + 4mf 2 ). 

As argued in [12] . for m 2 > ^ — \- 1, only the term with A+ is normalizable and 
only ^+ can be interpreted as the expectation value of an operator in a dual theory 
with dimension A+ . For = 4- + 1 > m 2 > ^- , both terms are normalizabla 11 !. However 
only one condensate can be turned on at a time in order to avoid instability in the 
asymptotic region of AdS. In this work we will confine ourselves to the case where only 
the mode with dimension A + is tuned on. 

4.1 Qualitative Discussion 

The hair-less solution to equations eqn (|4.21|) is given by, 

A t = Jl-r d - 2 ), Vi = 0, ^ 2 = (4.24) 



This is the normal phase solution that exists for all temperature. This solution is 
just the non-backreacted version of RN solution where back reaction of the gauge fields 
on the geometry is ignored. 

Next, let us consider the case when value of one of the scalars is set to zero. The 
problem then reduced to the well-studied holographic superconductors [13]. Let us 
assume that with the charge of the scalar fixed to unity, an instability in the scalar 
E.O.M. occurs at a critical value of fi = /ij 12 l. and a zero mode of the scalar field 
forms. If n is increased further, the zero mode condenses and a new phase with a 
non-trivial profile for ip is formed. This is the well known mechanism for holographic 
superconductors [13H15J . fi c is a increasing function of the mass square (m 2 ) of the 
scalar. As we vary the charge (e) of the scalars, p c is scaled to — . From the above 
discussion we may write down three solutions of the eqns (|4.2ip , 

• The hairless solution (eqn (]4.24|) ) exists for all value of p. And becomes locally 
unstable for fi > fi C 2. 

• The solution with ipx = and ip2 condenses. Exists for /i > /i C 2. 

• The solution with ip2 = and ipi condenses. Exists for fi > \i c \. 

Depending on the values of ei and e\, the second or third solution may become 
locally unstable. There may be new phases where both of the scalars condense, or it 
may also happen that one of the solutions (either the second and third solution) simply 
dominates over the other. 

To understand these issues in more detail we rewrite scalar EOM in ([4.2ip with a 



11 Note that the conformal mass is the case m 2 = —2, for d = 3. 

12 In a conformal theory, only the ratio fi/T matters. Increasing the chemical potential is equivalent to 
lowering the temperature. In the following, we choose to vary the chemical potential keeping the temperature 
fixed 
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generic background A t as, 

f / p 2 A 2 m 2 

r + Y + (Ij?- - y )</> = (4.25) 

I ^r 2 = ^ 

/(/VO -K//(r)^ = 0, \ v ^ f ^_ t2 ( (d-l)(d-3)r=» , r3(d-l)/' , e^| 



^e//(r) = -r - V "~ 1A r ; + -^TT^ + 



,2 



2/ T T / 

ci 2 ~ dr 

=*• ^2^ - K//(y)^ = 0, with dy = -^j. 

Here y — >• oo as r — )■ 1 and y — > as r — >• 0. In terms of this new variable y, the EOM 
of ?/> is rephrased as a potential problem on a semi infinite line (y : [0, oo)). Depending 
on the nature of the potential V e ff, there may exist a bound state of ifi. Such a bound 
state signifies an instability, suggesting that there may be new phases with non-trivial 
ip. In the absence of A t , the scalar ip may condense if and only if m 2 < tt? bf . This is the 
well known Bretinlohner-Freedman (BF) bound. Generically by choosing a sufficiently 
negative V e ff one may force ip to condense, e.g. if we use the normal phase value of 
At (eqn. I4.24p . and increase \i from zero, then eventually there will be zero mode of 
ip. Increasing fj, results in a bound state of ip signifying superconducting instability. 

Due to the Higgs mechanism the condensation of one scalar depletes the gauge field 
from the IR region of AdS space. This in turn reduces the amount of negative potential 
contributing to the condensation of another scalar. Hence we find that generically 
condensation of one scalar hinders the condensation of another^ However whether 
or not a condensation of one scalar field will completely stop the condensation of the 
other is a question of detail. 

From the above discussion it is natural to guess that, depending on whether fi C 2/^2 3> 
A*ci/ e i) the condensation of second scalar field dominates over the first one, or vice 
versa. Fortunately one may have some analytic arguments to understand these facts 
better and narrow down the range of possibilities. 

4.2 Analytic argument: 

Let us start by stating a simple lemma: if we define two potentials V\ and V% over 
the same domain and V\ > Vi then the lowest eigenvalue of V\ would be strictly 
greater than lowest eigenvalue of V 2 . This simple fact may be proven using variational 
argument, as follows. 

If the lowest eigenfunction in potential V\ is tpi with an eigenvalue Ai, 



Ai = J (-Vi< + vm 

> [{-ikl4 + V 2 il%) = \2 (4.26) 



The above derivation implies if V 2 does not have a bound state solution, then V\ 
cannot have one either. Also if the lowest eigenvalue mode for V 2 is a zero mode then 
Vi can not have a bound state or a zero mode. 
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This was observed in the case of non-Abelian holographic superconductor in [16] 
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Being in AdS complicates the picture but essence of the argument still survives. 
As we have discussed the scalar part of the Lagrangian may be written as a potential 
problem, 

Scalar = " J dr -1_ f m ^ - I^V + f\^f) 



dyU-^ + Ve ff (y)^) (4.27) 



dy 2 

Let us now get back to the two scalar problem. At first we will look at the case 
e 2 > e\. 



e 2 > e x case 

In this case we always have Vie// > V^//I 14 l. No matter which gauge field configuration 
we choose, this fact remains true. The above argument implies then that a zero mode 
of ipi may not form before a zero mode of "02 does, as /x is increased from a small 
normal phase value. Consequently ip2 condenses before ip\. 

Now the question is whether a zero mode of \j)\ may form in the phase where ?/> 2 
has already condensed. The answer is no, and the argument is as follows. It is to 
be noted that the 1^2 condensate is a zero mode of ip2- This is just the retelling of 
the fact that a non-trivial node-less solution exists for E.O.M. of ip2- Existence of 
such a solution comes primarily from numerics. The condensation of i/> 2 changes the 
gauge field profile and correspondingly change the potential problem. In this modified 
potential problem the condensate is a zero mode of ifa. This enables us to use the 
lemma proven at the beginning of this section. Since with the modified (i.e. modified 
by the condensation of ■0 2 ) gauge field V± e ff > V2 e ff and ip2 has a zero mode, it follows 
from the lemma that no zero mode or bound state of ipi exists. Hence in the phase 
with a tp2 condensation, tpl can not condense. This is in accordance to our general 
expectation: once ip2 condensed, it depletes the gauge potential, thus changing the 
effective potential for ip\ from one that allows a bound state, to one that does not. 

The slight modification of the above argument actually prevents any state with 
simultaneous condensation of ipi and V>2> as both may not have zero modes in any gauge 
field configuration. Moreover any phase where ip\ has a condensed will necessarily have 
a instability for ip2- Hence for e 2 > e%, any phase with non-zero if)% is dominant. The 
phase structure of the system is therefore the same as that of a single scalar holographic 
superconductor with the scalar field ip2- 

Next, let us consider the converse situation, where e 2 < e.\. 

e 2 < e x case 

The situation is more complicated here. Looking at the Lagrangian, one may naively 
expect that as we increase the chemical potential the scalar with more charge will al- 
ways eventually dominate. However the situation is more complicated as the potential 
V e ff in eqn (J4.25P diverges like -3- near the boundary y = 0. In this case the above ar- 
gument is essentially reversed and the mass dependent potential part possibly becomes 
important as we turn on more chemical potential. This concurs with the numerical 
finding discussed next. 

14 Remember ml > m\. 
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4.3 Numerical Analysis 

Here we study the eqn (|4.2ip numerically. We confine ourselves to the case with 
m\ = and m\ = rn1 on f orma i = — 2 (see p2]). With just one of the scalar turned 
on and the charge of it set to unity we know that the scalar condenses as the value 
of [i is increased from zero. Let us assume that without any interface from another 
scalar, scalars condense at // = \x c . If we set the charge of the scalars to unity then 
/i cl m= 7.8 « 1.9 / u C 2(w 4.0). As the charge (e) of the scalar varies from unity, \i c is 
scaled to /i c /e. 

As discussed before in the regime &2 > &i the phase diagram is same as that of a 
model with only second scalar. The same picture holds in the regime e\ < jf^£2 ~ 1.9e2 
and the second scalar dominates the picture and there is no condensation of the first 
scalar. 

As e2 decreases beyond —e\ the first scalar condenses before the second scalar. As 

fJ'C 1 

we increase the boundary value of chemical potential more the second scalar condenses 
at a chemical potential [i = [i c %. The condensation of second scalar increases with 
chemical potential. As argued in the previous section increasing chemical potential 
possibly makes the mass dependent potential term more important and hence the 
second scalar tends to dominate the picture. As \x is increased further the condensation 
of first scalar decreases and goes to zero at /J, = j2 c i ■ Although the first scalar condenses 
before the second one, one may consider the solution where the expectation value of 
first scalar is kept to zero and let the second scalar condenses. As chemical potential is 
increased, there will be some value jx c \ when the first scalar will develop a zero mode. 
Numerically fi C 2 is close to \x c \ and consequently one phase quickly changes to another. 
In this case there are three apparently second order transition at the increasing order 
of fi (see fig. Oil: 

1. A second order phase transition at // = ix c \ when first scalar condenses. 

2. A second order phase transition at /x = /2 C 2 when second scalar start condensing. 
The resulting phase has both scalar condensates. 

3. A when the condensation of first scalar goes to zero at [i = fj, c i. 

Numerically it seems that as e2 is decreased even further such that e^ < 1.15^-ei ~ 
2.2ei the second scalar never condenses and the resulting phase diagram is same as 
that of model with only first scalar. Whether this is a numerical artefact or not needs 
to be investigated in more detail. 

In principle there may exist first order transitions in our system. The particular 
case we have discussed seems to have none. A detailed study with other possible values 
of mass and dimension is left out for a future work. 

4.4 Effect of Scalar Interactions 

Similar to Lifshitz case discussed above, turning on a direct interaction between scalars 
should help scalars to condense simultaneously, in case that interaction is repulsive. It 
is then expected that the regime of co-existence of phases will be enhanced. We expect 
an opposite effect for a repulsive interaction. Detailed study of these issues are left for 
a future work. 
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Figure 5: Plot of condensate with two scalars with e\ = 1.95 and e<i = 1. The curve on the 
right is for the second scalar and it starts from [i = £l c2 . The left one is first scalar it starts 
from fi = fii and ends at \i = Ji c \. In the range \x C 2 < fj, < jl c i both of the scalar condenses. 
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